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ABSTRACT. There are two arylliydrocarbon receptor (AliR) isoforins in birds, AliRl and AhR2. The varying sensitivity of AhR is reported 
to be related to two critical amino acids at positions 325 and 381 in the AhRl ligand-binding domain. In this study, seven avian species 
whose in vivo dioxin sensitivity was known, and 1 3 species with no data regarding their in vivo dioxin sensitivity were examined. The two 
critical amino acids in the ligand-binding domain were investigated in avian species, and the results were compared with the taxonomy or 
phylogenetic trees for the bird AhR proteins. We found that the two critical amino acids did not correlate with the taxonomy or phylogeny 
of these proteins, suggesting that dioxin sensitivity was independent of taxonomy. 
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Environmental pollutants, such as 2,3,7,8-tetrachloro- 
dibenzo-/7-dioxin (TCDD), halogenated aromatic hydro- 
carbons (HAHs) and polycyclic aromatic hydrocarbons 
(PAHs), can induce serious toxicity in avian species. The 
types of toxicity induced include teratogenic, immunotoxic 
and reproductive toxicity [1, 2, 19, 23, 31]. The dioxin con- 
centration required to induce these types of toxicity varies 
significantly between birds [3]. The large difference in 
dioxin sensitivity among avian species is reported to be de- 
pendent on the arylhydrocarbon receptor (AhR) protein [11, 
17], which has a role in the induction of toxicity [10, 21, 27]. 

AhR is a basic-helix-loop-helix/Per Arnt Sim (PAS) fam- 
ily protein and a transcription factor activated by ligand 
binding [6]. When not bound to a ligand, AhR remains in 
the cytosol, fonning a complex with heat shock protein 90 
(HSP90), AhR-associated protein (XAP2 or ARA9) and p23 
[7, 24]. Once bound to a ligand, AhR is translocated to the 
nucleus [32] where it forms a heterodimer with an AhR nu- 
clear translocator (Arnt), which then binds to the xenobiotic 
responsive element (XRE) [20, 28]. After binding to XRE, 
transcription of the CYPlAl, CYP1A2 and AhR repressor 
(AhRR) genes is activated [13]. 

Avian species have two AhR isoforms, AhRl and AhR2 
[33, 34], whereas most mammals possess only one. The 
dominant isoform of AhR differs among bird species [18], 
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and there are large differences in fiinction, even within the 
same AhR isoform. For example, although avian AhRl s are 
highly conserved (>90%) among species, there are large 
interspecies differences in their sensitivity to dioxins, which 
can be explained by differences in their ligand-binding af- 
finities and transactivation abilities. 

It is reported that AhR sensitivity can be predicted from 
the two amino acids at positions 325 and 381 of AhRl [17]. 
Chicken is well known to be the only avian species which 
has a sensitive type of AhR, however, the sensitivity of AhRs 
in broad avian species is still unclear. In this study, several 
kinds of avian species, which were chosen from phyloge- 
netic tree of bird, were investigated to determine their dioxin 
sensitivity. The amino acid sequences of AhRl and AhR2 
were determined for each species, and compared to their 
taxonomic and phylogenetic classifications. 

MATERIALS AND METHODS 

Animals: Bird species analyzed in this study were se- 
lected considering clade of phylogenetic tree based on DNA 
sequences [12]. A one-year-old female blue-eared pheasant 
{Crossoptilon auritum), a male ruddy shelduck {Tadorna 
ferrugined), a one-year-old male mallard {Anas platyrhyn- 
chos), two male great horned owls {Bubo virginianus), two 
male and one female bar-headed goose {Anser indicus), one 
male Indian peafowl {Pavo cristatus), one 12-year-old male 
goose {Anser anser), one 19-year-old female black-headed 
ibis {Threskiornis melanocephalus), one female swan goose 
{Anser cygnoides), two male snowy owls {Bubo scandiacus), 
one female Chilean flamingo {Phoenicopterus chilensis), one 
eight-year-old male Humboldt penguin {Spheniscus hum- 
boldti), one female cape barren goose {Cereopsis novaehol- 
landiae) and one gender-undetermined black-crowned night 
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Table 1 . Primers for avian AhR 




Forward 5'-sequence-3' 


Reverse 5'-sequence-3' 


Avian AHR Full 
Avian AHRl-1 
Avian AHR 1-2 
Avian AHR 1-3 
Avian AHR 1-4 


CAGGATGAACCCCAATGTCAC 

GGATGAACCCCAATGTCACCTA 

TCATCTGCAGGTTACGATGCCT 

TGCCCTTCATGTTTGCCACTGGTGA 

CAGCTCTGTCAAAAGATGAAA 


GTCACATAAATCCACTAGATGCCAAA 

ATCGTCCTTGAAAATTCATA 

ACACAGACTCATCTTGCCTTA 

TCCAATTTGTGAACATCCCAT 

TTACATAAATCCACTAGA 



heron {Nycticorax nycticomx), were provided by Maruyama 
Zoo (Sapporo, Japan). These animals died due to accidental 
injury or disease, such as enteritis. Their livers were immedi- 
ately frozen in liquid nitrogen and stored at -80°C until use. 
All experiments using animals were performed according 
to the guidelines of the Hokkaido University Institutional 
Animal Care and Use Committee. 

cDNA cloning and sequencing of AhR: Hepatic total 
RNA was isolated using TRI reagent (Sigma-Aldrich, St. 
Louis, MO, U.S.A.) and reverse -transcribed to cDNA using 
Oligo(dT). Partial AhRl and AhR2 DNA sequences were 
amplified by PCR using the primers listed in Table 1 . The 
PCR parameters for the amplification of AhRl were as fol- 
lows: 94°C for 2 min, then 94°C for 30 sec, 66°C for 45 sec 
and 72°C for 3 min for 35 cycles, followed by 72°C for 5 
min. PCR products were subject to direct-sequencing using 
primers Avian AhRl- 1 to 4 and an annealing temperature of 
50°C; the forward and reverse AhR2 primers were used in 
pairs with an annealing temperature of 63°C. For obtaining 
DNA sequence information for the ligand-binding domain, 
the primers: AhR2-LBD-F 5'-TCTCCAGACAAAGCA- 
CAAGCTGGAC-3' and AhR2-LBD-R 5'-GTACAG- 
GACTGCTTCCCCCGTG-3' were used. Reproducibility of 
sequence was confirmed at least three times. 

Phylogenetic tree and amino acid sequence alignment 
of AhRl and AhR2: DNA sequences of avian AhRl were 
aligned by CLUSTAL W using Molecular Evolutionary 
Genetics Analysis (MEGA) 5 [30]. The accession numbers 
for the sequences included in this analysis were: ostrich 
AhRl (AB820092), blue eared-pheasant AhRl (AB820094, 
AB820095), Indian peafowl AhRl (AB820097), swan goose 
AhRl (AB820099), bar-headed goose AhRl (AB820101), 
goose AhRl (AB820103), mallard AhRl (AB820105, 
AB820106), ruddy shelduck AhRl (AB820108), cape bar- 
ren goose AhRl (AB820110), black-headed ibis AhRl 
(AB820112), Humboldt penguin AhRl (AB820113), Chil- 
ean flamingo AhR 1 ( AB 820115, AB 820116), black-crowned 
night heron AhRl (AB820118, AB820119), snowy owl 
AhRl (AB820121), great horned owl AhRl (AB820122, 
AB820123), peregrine falcon AhRl (AB560859), common 
cormorant AhRl (AB 109545), black-footed albatross AhRl 
(AB106109), common tern AhRl (AF192503), chicken AhR 
(AF192502), ostrich AhR2 (AB920093), blue eared-pheas- 
ant AhR2 (AB820096), Indian peafowl AhR2 (AB820098), 
swan goose AhR2 (AB820100), bar-headed goose AhR2 
(AB820102), goose AhR2 (AB820104), mallard AhR2 
(AB820107), ruddy shelduck AhR2 (AB820109), cape 
barren goose AhR2 (AB820111), Humboldt penguin AhR2 



(AB820114), Chilean flamingo AhR2 (AB820117), black- 
crowned night heron AhR2 (AB820120), great horned owl 
AhR2 (AB820124), peregrine falcon AhR2 (AB560860), 
black-footed albatross AhR2 (AB106110), common cormo- 
rant AhR2 (AB287294) and chicken AhR2 (XM421887). 
Human AhR (LI 9872) was added as an outgroup. For AhR2, 
an alignment was performed using -250 bases of DNA of the 
ligand-binding domain and excluded areas containing gaps. 
A phylogenetic tree was constructed by the Maximum likeli- 
hood method based on MEGA5 program. The bootstrap con- 
sensus tree inferred from 500 replicates is taken to represent 
the evolutionary history of the taxa analyzed [9]. Tamura- 
Nei model [29] was applied into the nucleotide sequences. 

RESULTS 

Two critical amino acids in the AhRl and AhRl proteins: 
Based on the two critical amino acids in the ligand-binding 
domain of AhRl [17], the avian species we examined could 
be divided into three groups. The first group, with amino ac- 
ids 325-Ile and 381-Ser, consisted of the ostrich and chicken. 
The chicken is reported to be a highly sensitive species to 
TCDD [14, 17, 18, 34]. The blue -eared pheasant, Indian 
peafowl, black-footed albatross and swan goose composed 
the second group, possessing amino acids 325-Ile and 381- 
Ala. Other avian species, including the bar-headed goose, 
goose, mallard, ruddy shelduck, cape barren goose, snowy 
owl, great horned owl, peregrine falcon, black-headed ibis, 
Humboldt penguin, Chilean flamingo, common cormorant 
and the black-crowned night heron, belonged to the last 
group, harboring amino acids 325-Val and 381 -Ala (Figs. 1 
and 2; Table 2). 

The species could also be divided into three groups ac- 
cording to the amino acids in the AhR2 ligand-binding 
domain. The group possessing the amino acids 325-Leu 
and 381-Ala comprised the ostrich, blue-eared pheasant and 
Indian peafowl. The chicken is the only one species to con- 
stitute the group of species to possess 325-Val and 381-Ser 
AhR2. Other avian species belonged to the group harboring 
the amino acids 325-Val and 381-Ala (Figs. 1 and 2; Table 
2). 

AhRl andAhRl phylogenetic analyses: The phylogenetic 
tree constructed from the avian AhRl sequences indicated 
that the ostrich was distinct from the other avian species. 
Also, Galliformes, including chicken, pheasant and peafowl, 
Ciconiiformes, including albatross, tern, cormorant, pen- 
guin, flamingo, falcon and heron and Strigiformes, including 
the snowy owl and great horned owl, grouped close together 
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Fig. 3. Phylogenetic analysis of avian AhRl. DNA sequences of AliRls were aligned by CLUSTAL W using the 
MEGA5 program. Human AhR (LI 9872) was added as an outgroup. Alignment was performed with a length of 
about 2,000 bases, including the functional domains, PAS-A, PAS-B and the Q-rich domains and excluded the 
regions containing gaps. The phylogenetic tree was constructed by ML method using MEGA5. The number of 
bootstrap replications was set to 500. The percentages of replicate trees in which the associated taxa clustered 
together in the bootstrap test (500 replicates) are shown next to the branches. Tamura-Nei model was applied in 
nucleotide sequences. Swan goose, bar-headed goose, goose, ruddy shelduck, cape barren goose and black-headed 
ibis are not shown in this figure. 



Table 2. The critical amino acids in ligand-binding domains of AhRl and AhR2 with in vivo sensitivity 





325 


AhRl 
381 


325 


AhR2 

381 


Order 


In vivo Sensitivity 


Ostrich 




S 


L 


A 


St 




Chicken 




S 


V 


S 


G 


High [16] 


*Blue-eared Pheasant 




A 


L 


A 


G 


Middle [3, 22]'-'» 


*Indian Peafowl 




A 


L 


A 


G 




*Swan Goose 




A 


V 


A 


A 




*Bar-headed Goose 


V 


A 


V 


A 


A 




* Goose 


V 


A 


V 


A 


A 


low [5] 


♦Mallard 


V 


A 


V 


A 


A 


low [3, 5, 16] 


*Ruddy Shelduck 


V 


A 


V 


A 


A 




*Cape Barren Goose 


V 


A 


V 


A 


A 




Black-footed Albatross 


1 


A 


V 


A 


C 




Common Cormorant 


V 


A 


V 


A 


C 


low [26]'" 


Peregrine Falcon 


V 


A 


V 


A 


c 


low [15]" 


Common Tern 


V 


A 






c 


low [4, 15] 


*Black-headed Ibis 


V 


A 


V 


A 


c 




*Humboldt Penguin 


V 


A 


V 


A 


c 




*Chilean Flamingo 


V 


A 


V 


A 


c 




*Black-crowned Night Heron 


V 


A 


V 


A 


c 




* Snowy Owl 


V 


A 


V 


A 


Sg 




*Great Homed Owl 


V 


A 


V 


A 


Sg 




The two amino acids in AhRl and AhR2 at positions 325 and 38 1 of each avian species were indicated. I: 
isoleucine, S: serine. A: alanine, V: valine, L: leucine. * indicates the species whose AhRs are cloned and 
sequenced in this study. Abbreviations for each order are; St: Stnithioniformes, G: Galliformes, A: Anseri- 



formes, C: Ciconiifonnes and Sg: Strigiformes. a) In vivo sensitivity of common pheasant or ring-necked 
pheasant (Phasiamis colclncus). b) vivo sensitivity of double-crested cormorant (Phalacrocorax auritus). 
c) vivo sensitivity of American kestrel (Falco spanerius). 
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Fig. 4. Phylogenetic analysis of avian AhR2. DNA sequences of AhR2s were aligned using the by CLUSTAL 
W using MEGA5 program. Human AhR (L19872) was added as an outgroup. Alignment was performed at 
a length of about 200 bases, a region of the ligand-binding domain, and also excluded regions containing 
gaps. The bootstrap consensus tree inferred from 500 replicates is taken to represent the evolutionary history 
of the taxa analyzed. The percentages of replicate trees in which the associated taxa clustered together in the 
bootstrap test (500 replicates) are shown next to the branches. The phylogenetic tree was constructed by ML 
method using MEGA5. Tamura-Nei model was applied in nucleotide sequences. Snowy owl and black-headed 
ibis are not shown in this figure. 



on the tree. These groupings were consistent with the taxo- 
nomic groupings of the avian species (Fig. 3). 

In the case of AhR2, the phylogenetic tree was also Hnked 
to the taxonomy of the species with the exception of the 
ostrich and great horned owl. The orders Galliformes, Ci- 
coniiformes and Anserifonnes in phylogenetic tree were all 
assembled the same as that of AhRl (Fig. 4). 

DISCUSSION 

In mammals, several factors have been reported to de- 
termine the ligand-binding affinity or dioxin sensitivity to 
AhR-induced toxicity. In C57BL/6 and DBA/2 mice, large 
differences in AhR ligand-binding affinity result from AhR 
point mutations at codon 375 in the ligand-binding domain 
[8]. Similarly, Han/Wistar rats, which are insensitive to di- 
oxin, harbor a point mutation at position 497 in AhR [25]. 

Avian species are distinct in that they possess two AhR 
isoforms, AhRI and AhR2, and the type of AhR dominantly 
expressed varies among avian species [18]. Most species 
dominantly express AhRI, and AhRI ligand-binding affinity 
is reported to directly correlate with CYPIA transactivation 
ability [18]. Critical mutations which decide dioxin sensitiv- 
ity have been found in avian species at positions 325 and 381 
of AhR I [17]. Head et al. [ 1 4] showed these two amino acids 
are effective for predicting avian dioxin sensitivity. In the 
study, avian species are divided into three groups according 
to the key amino acids in AhRI . The most sensitive group is 



with AhRI of 325-IIe and 381-Ser, middle for 325-IIe and 
381-Ala, and the least sensitive for 325-Val and 381-Ala. In 
this study, the blue-eared pheasant, goose and mallard are 
newly classified according to the two amino acids, and we 
found that the grouping corresponds with in vivo sensitivity 
also in these cases. 

In the current study, we investigated whether similar 
avian species would have similar levels of dioxin sensitivity. 
Indeed, the phylogenetic trees constructed from the amino 
acids sequences of AhRs gave results in-keeping with the 
evolutionary history of these proteins [12]. Unexpectedly, 
even though the amino acids sequences of AhRs highly re- 
flected the taxonomy, the identity of the two critical amino 
acids, at positions 325 and 381, did not correspond to the 
phylogenetic trees of AhR or taxonomy. In fact, these two 
amino acids were not conserved in the orders, Galliformes 
and Ciconiiformes [14, 34]. Therefore, our new findings 
suggested that these key amino acids at positions 325 and 
381 are independent from the other amino acids sequences 
of AhRs, so that they cannot be predicted from the phylo- 
genetic tree or from taxonomy. That is, the ligand-binding 
affinity or the dioxin sensitivity of each avian AhR protein 
cannot be determined from the taxonomy. 

The amino acid sequence of the AhRI ligand-binding do- 
main was highly conserved among different species. In the 
ostrich, multiple amino acid changes were found throughout 
the ligand-binding domain, compared with other avian spe- 
cies. This ostrich AhRI was reported to possess high trans- 
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activation ability in our previous study, similar to chicken 
AhRl [11]. 

Regarding the two critical amino acids, avian species 
harboring isoleucine at position 325 in AhRl were chicken, 
peafowl, pheasant, albatross, swan goose and ostrich. All 
species in the order Galliformes examined in this study, in- 
cluding chicken, peafowl and pheasant, possessed isoleucine 
at position 325. However, albatross and swan goose were the 
only species to possess this amino acid in their respective 
orders, Ciconiiformes and Anseriformes. Only two of the 
avian species, chicken and ostrich from the orders Gallifor- 
mes and Struthioniformes, respectively, possessed a serine at 
position 381. Species within the order Galliformes, such as 
peafowl and pheasant, did not harbor this amino acid. Taken 
together, these findings indicate that the two critical amino 
acids are independent of taxonomy or even phylogeny of 
the full-length amino acid sequence of AhRl . Therefore, we 
conclude that it is difficult to predict the dioxin sensitivity of 
avian species from taxonomy or evolutionary history. 

In the case of avian AhR2, the amino acid sequence of 
the ligand-binding domain was not as highly conserved 
as that of AhRl. In terms of the two critical amino acids 
in the ligand-binding domain, the only species to possess 
381-Ser in AhR2 was the chicken. This amino acid was 
not conserved in the AhRl and AhR2 of ostrich. In addi- 
tion, we identified the amino acid leucine at position 325 
in pheasant, peafowl and ostrich. This amino acid has not 
previously been reported at this position, and its correspond- 
ing AhR fiinction is therefore unknown. It will be of interest 
to investigate the fiinction or ligand-binding affinity of this 
type of AhR protein. Future researches are also required to 
fully investigate the avian AhR2 protein and its role in avian 
dioxin sensitivity. 

In conclusion, the two critical amino acids at positions 
325 and 381 in the ligand-binding domain of AhR were 
investigated in several bird species, and the results were 
compared with the taxonomy or phylogenetic trees for the 
AhR proteins. The two critical amino acids did not correlate 
with the taxonomy or phylogeny of these proteins, and di- 
oxin sensitivity was independent of taxonomy. 
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